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I.  RESEARCH  PROGRAM 

Major  accompli sh merits  for  the  report  period  are  described  in  Section  II 
below.  Reports  detailing  these  accomplishments  will  be  included  in  the  final 
report  for  this  contract. 

II.  MAJOR  ACCOMPLISHMENTS 

A.  Electronic  f-numbers 

Programs  for  the  calculation  of  electronic  f- numbers  have  been  written  in 

two  general  parts.  M.  Yosh inline  has  completed  and  debugaed  a  program  for 

the  calculation  of  matrix  elements  of  r  between  Cl  wave  functions  for  different 

electronic  states.  The  program  permits  the  KO's  used  for  the  construction  of 

the  final  state  wave  function  to  be  different  from  those  used  for  the  initial 

state  wave  function.  The  Cl  expansions  involved  may  be  quite  long;  for 
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example,  for  the  X  It  and  A  X  states  of  Li 0  approximately  1500  symmetrized 
configurations  and  -  4500  determinants  were  used  for  each  state.  Because 
of  these  long  expansions,  the  calculation  of  ^ jx.r |¥p>  may  become  reasonably 
time  consuming.  To  reduce  the  computational  time,  thresholds  for  the  selection 
of  only  a  part  of  the  total  number  of  configurations  and  a  part  of  the  total 
number  of  matrix  elements  between  configurations  have  been  included  as  an  essen¬ 
tial  part  of  the  program,  for  the  second  general  part  of  this  work,  B.  Liu 
has  extended  his  vibration-rotational  analysis  program  to  calculate  electronic 
band  f-numbers  and  Frank-Condon  factors.  This  program  uses  as  input  the  electronic 
transition  matrix  elements  of  r  generated  by  the  first  program. 

M.  Yoshimine  has,  using  the  programs  de  cribed  above,  obtained  electronic 
band  f-numbers  for  the  A^£+  transition  of  Li 0  end  the  X'V*-* A*HI  transition 
of  A10. 

B.  Sigma  States  of  A10 

B.  Liu  anc  P.  S.  Bagus  attempted  to  use  the  multi-configuration  SCF  (MCSCF) 
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method  to  stucy  the  X  T,  state  of  A10.  This  state  has  an  anomolous  behavior 
in  the  one-corfiguration  SCF  approximation  vhich  caused  considerable  difficul¬ 
ties  in  the  selection  of  a  reference  state  for  subsequent  Cl  calculations.  It 
was  hoped  thal  the  MCSCF  method  would  prov'de  an  unambiguous  reference  slate 
avoiding  these  difficulties.  Unfortunately,  the  self-consistent  field  itera¬ 
tions  for  the  MCSCF  orbitals  failed  to  converge. 
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In  order  to  gain  insight  .into  the  B  T 
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performed  limited  Cl  calculations  on  'l  states.  These  calculations  involved, 
beyond  the  orbitals  occupied  in  the  y?'7.  SCI'  calculations,  two  additional 
sinma  and  one  additional  pi  orbital.  The  total  number  of  configurations  was 
140  and  both  the  lowest  ( X^>:  )  and  second  lowest  (B  l  )  roots  were  obtained, 
i he  energy  separation  of  the  two  stales  is  10  percent  larger  than  that  observed 
experimentally.  Electronic  f-nutnbers  will  be  calculated  with  these  wave  func¬ 
tions  and  compared  to  those  obtained  by  Michels  [J.C.P.,  5£,  665  (1972)]. 


C.  Electronic  States  of  FeO 

P.  S.  Bagus  has  completed  atomic  correlation  energy  calculations  on  Fe 
and  0  using  the  same  set  of  excitations  as  used  for  the  V  state  of  FeO. 

H* 

These  calculations  give  a  dissociation  energy  of  1.2eV  for  the  >:  state. 
However,  because  of  the  limited  number  of  molecular  orbitals  used  iri  %>7  Cl 
calculation,  the  treatment  of  the  separated  atoms  is  expected  to  be  far  better 
than  for  the  molecule.  To  avoid  this  problem,  pair  correlation  energy  calcu¬ 
lations  on  1  he  V  state  using  essentially  a  full  basis  set  of  HO 1 s  were  per- 
thr«p  pninfs  around  the  computed  minimum  of  the  potential  curve. 
These  calculations  indicate  that  the  dissociation  energy  is  below  2eV.  liris 
is  much  lower  than  the  previous  estimate  but  still  short  of  the  ~4eV  observed 
for  the  ground  state. 

Earlier  calculations  on  the  singly  ionic  V  state  have  been  extended. 

SCF  calculations  have  been  performed  for  7>V  states  on  models  of  neutral 
and  of  doubly  ionic  systems.  In  addition,  limited  Cl  calculations  have  L een 
performed  which  mix  all  of  these  types  of  configurations.  The  doubly  ioiic 
configuration,  iri  particular,  was  found  to  have  an  important  effect  on  tie 
behavior  of  this  state.  In  addition  SCF  ca  culations  on  a  singly  ionic  model 
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for  the  l~  state  ‘were  attempted. 

For  most  states,  SCF  convergence  problems  were  encountered  and  for  several 
states  convergence  could  not  be  obtained,  ihe  use  of  adjusting  equivalent, 
in  principle,  forms  of  Fock  operators  to  obtain  convergence  was  investigated 
with  limited  success. 


D.  Polyatomic  Matrix  Elements 

A.  D.  McLean  has  worked  on  the  analysis  of  the  evaluation  of  matrix 
elements  between  Slater  type  basis  functiors  for  molecules  of  general 
geometry.  A  t eport  on  methods  for  evaluation  is  near  completion. 


Program  Use  and  Large  Scale  Cl  Methodology 


s/ 


E.  ALCHEMY  Documentation: 

P.  S.  Banns  has  completed  an  article  describing  the  theory  of  open 
shell  energy  expressions  and  their  use  in  SCI  calculations.  A.  D.  McLean 
continues  work  on  an  ALCHEMY  "User  Manual"  including  detailed  reports  on 
the  methods  currently  available  for  use  with  this  program  system. 

A.  D.  McLean  and  B.  Liu  have  developed  a  scheme  for  the  classification 
of  configuration  state  functio,.  forming  the  11-particle  basis  of  Cl  calcu¬ 
lations.  The  classification  includes  the  order  of  perturbation  theory 

in  which  a  configuration  at.c  function  makes  a  correction  to  a  given 
zeroth  order  function.  They  have  shown  how  to  minimize  the  number  of  con¬ 
figuration  state  functions  having  a  non-zero  matrix  element  with  members 
of  the  zeroth  order  function.  This  work  is  of  importance  in  reducing  the 
size  of  Cl  calculations  designed  to  yield  Lome  prescribed  accuracy,  at 
the  same  time  it  makes  available  a  language  for  interpreting  a  wave  function 
and  is  cf  assistance  in  preparing  Cl  calculations.  A  report  on  this  work  is 
in  preparation. 

F.  Methods  for  Large  Scale  Cl  Calculations 

B.  Liu  has  worked  on  methods  for  applying  Rayleigh-Schrddinger  perturbation 
theory  as  a  guide  for  selecting  configurations  in  large  Cl  calculations.  He 
has  also  used  |erturbation  theory  as  a  method  to  obtain  approximate  natural 
orbitals.  The  method  was  studied  and  applied  in  calculations  on  the  low-lying 
states  of  CH.  Effective  methods  for  selection  of  configurations  and  approxi¬ 
mate  natural  orbitals  are  essential  to  make  accurate  Cl  calculations  tractable. 

III.  Problems  Encountered 
None 

IV.  Fiscal  Sta  :us 

Research  Libor  and  Research  Burden  expenditures  under  this  contract  ior 
the  period  June  1C,  1971,  to  May  19,  1972,  were  $68,178;  computer  time  aid 
support  time  '.Larges  for  the  same  period  we 'e  $13,112. 
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Action  Required  by 
None. 


P.  S.  Bagus 

Principal  Investigator 
PSBrjni 
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tne  Government 


